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DEVICE AND METHOD FOR ESTIMATING PHYSICAL QUANTITY AND DEVICE AND 

METHOD FOR ABS CONTROL 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to a physical quantity estimation device for 
estimating a physical quantity relating to the ease of slippage between a wheel and the 
surface of a road, and particularly to a physical quantity estimation device capable of 
accurately estimating a physical quantity relating to the ease of slippage between a 
road surface and a wheel even when an unsprung resonance phenomenon is 
generated, and to an ABS control device for controlling braking pressure using this 
physical quantity estimation device. 
Description of the Related Art 

Japanese Patent Application Laid-Open (JP-A) No. 10-114263 discloses an anti- 
lock braking system (ABS) control device for performing, on the basis of a wheel 
deceleration dynamic model, a calculation to estimate the torque gradient, i.e. the road 
surface u gradient (the ratio between the road surface At and the slip) from the time 
series data of a wheel velocity signal, and for then controlling the manipulation variable 
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of the braking force so that the road surface m gradient is 0. The "road surface p. 
gradient" herein means a coefficient of friction u between the wheel and the road 
surface. 

When a vehicle travels from an icy road surface onto a dry road surface while the 
braking force is being controlled so as to be at the maximum amount by the above 
described ABS control device, then, as is shown in Rg. 1A, an oscillation phenomenon 
is generated in the wheel velocity. This wheel velocity oscillation phenomenon causes 
the accuracy of the u gradient estimation to be reduced. As is snown in Rg. IB, in 
some cases, a delay in the adaptation of the ABS control is also generated. 

The present inventors discovered that this phenomenon of oscillation in tiie 
wheel velocity is due to resonance movement between the tires and the suspension in 
the longitudinal direction of the vehicle. 

SUMMARY OF THE INVENTION 

The present invention is based on the above discovery and aims to provide a 
physical quantity estimating device capable of accurately estimating a physical 
quantity relating to the ease of slippage between a wheel and the surface of a road 
even when longitudinal resonance movement is generated between the tires and the 
suspension. 



In order to achieve the above aim, -the present invention comprises: wheel 
velocity detection means for detecting a wheel velocity and outputting a wheel velocity 
signal; and physical quantity estimation means for identifying a parameter of a 
physical model representing unsprung resonance characteristics from the wheel 
velocity signal, and estimating a physical quantity relating to ease of slippage between 
a road surface and a wheel from the identified parameter. 

According to the present invention, the wheel velocity is detected by the wheel 
velocity detection means and a wheel velocity signal is output. The physical quantity 
estimation means identifies the parameter of a physical model representing the 
unsprung resonance from the wheel velocity signal and estimates a physical quantity 
relating to the ease of slippage between the road surface and the wheel from the 
identified parameter. 

The physical model of the present invention can be a resonance model which has 
the unsprung resonance frequency as the undamped natural frequency. In this case, 
the physical quantity estimation means estimates a physical quantity relating to the 
ease of slippage between the road surface and the wheel on the basis of the damping 
ratio of the resonance model. 

In the present invention, because a parameter of a physical model representing 
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unsprung resonance characteristics is identified and a physical quantity relating to the 
ease of slippage between the road surface and a wheel is estimated from the identified 
parameter, the physical quantity can be accurately estimated even when an unsprung 
resonance phenomenon is generated. 

The present invention is capable of being used in an ABS control device for 
controlling braking pressure on the basis of the estimated physical quantity. . 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figs. 1A and IB show the results of an experiment using a conventional ABS 
control device when a vehicle travels from an icy road surface onto a dry road surface. 
Fig. ia is a line graph showing the relationship between the wheel velocity and time, 
while Fig. IB is a line graph showing the relationship between the braking wheel 
cylinder pressure and time. 

Fig. 2 is a structural diagram showing a suspension - tire resonance model. 

Figs. 3A and 3B show relationships between the road surface ti gradient and 
the suspension r tire resonance characteristics. Fig 3A is a line graph showing the 
relationship between the frequency and the gain, while Fig. 3B is a line graph showing 
the relationship between the frequency and the phase. 

Figs. 4A to 4G show the results of identifying the damping ratio in the experiment 



in which a vehicle with the ABS control device according to the present embodiment 
travels from an icy surface onto a dry surface. Fig. 4A is a line graph showing changes 
over time in the vehicle speed and wheel velocity, Fig. 4B is a line graph showing 
changes over time in the master cylinder pressure and wheel cylinder pressure, and 
Fig. 4C is a line graph showing changes over time in the damping ratio. 

Figs. 5A to SCshow the results of ABS control using damping ratio identification 
values according to the present embodiment Fig. 5A is a line graph showing changes 
over time in the vehicle speed and wheel velocity, Fig. 5B is a line graph showing 
changes Over time in the master cylinder pressure and wheel cylinder pressure, and 
Fig. 5C is a line graph showing changes over time in the damping ratio 

Rg. 6 is a block diagram in which the physical quantity estimation device 
according to the embodiment of the present invention is used in an ABS control device. 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 
A detailed description of an embodiment in which the physical quantity 
estimation device of the present invention is used in an ABS control device with 
reference to the Figures is given below. 

As is shown in Rg. 6, The ABS control device according to the present 
embodiment is provided with a wheel velocity sensor 10 which detects the wheel 



velocity of each of the four wheels of a vehicle arid outputs a signal such as an electric 
signal or a numerical signal corresponding to the wheel velocity as a wheel velocity 
signal. A bandpass filter is connected to the output terminal of the wheel velocity 
sensor 10. The bandpass filter removes signals of frequency bands which are 
unrelated to the unsprung resonance from the wheel velocity signal, and passes on 
only signals of frequency bands which are related to unsprung resonance. • 

An ju gradient estimating device 14 is connected to the output terminal of the. 
bandpass filter 12. The u gradient estimating device 14 uses an online identification 
method to identify, as a parameter, the damping ratio of a second order resonance 
model closely resembling the suspension - tire resonance model which is a physical 
model representing unsprung resonance characteristics from the signals output from 
the bandpass filter 12. The u gradient estimating device 14 then estimates the road 
surface p. gradient (referred to below simply as the n gradient), which is one of the 
physical quantities relating to the ease of slippage between the road surface and the 
wheel, from the identified damping ratio. Note that the wheel velocity sensor 10, the 
bandpass filter 12, and the u gradient estimating device 14 form the physical 
quantity estimation device of the present invention. 

Ah ABS control section 16 is connected to the output terminal of the # gradient 
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estimating device 14. The ABS control section 16 calculates the manipulation variable 
of the braking force acting on each wheel so. that the p. gradient is 0 or a small 
positive value near 0. An ABS control valve 18, which controls the wheel pressure, is 
connected to the output terminal of the ABS control section 16. 

According to the present embodiment, the wheel velocity signal output from the 
wheel velocity sensor 10 is input to the bandpass filter 12, and only signals related to 
the unsprung resonance are output from the bandpass filter 12. Signals in the 
frequency band related to the unsprung resonance are input to the il gradient 
estimating device 14 and the u gradient is estimated in the method explained below. 
Next, the manipulation variable of the braking force acting on each wheel is calculated 
by the ABS control section 16 so that the estimated U gradient is 0 or a small positive 
value near 0. The braking master pressure is then controlled by the ABS control valve 
18 being controlled in accordance with this manipulation variable. 

Next, the principle of the u gradient estimation of the present embodiment will 
be explained. Firstly, the relationship between resonance characteristics and ul 
gradient is explained. The resonance characteristics are those obtained when using a 
suspension - tire resonance model which takes into account the unsprung 
characteristics is postulated as a model for expressing the wheel velocity resonance 
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characteristics. As is shown in Rg. 2, this suspension - tire resonance modei is formed 
by connecting a tire T to a vehicle body B via a suspension S. 

The equation (the resonance model between the tires and the suspension in the 
longitudinal direction of the vehicle) of motion expressing disturbances around the 
equilibrium point in Fig. 2 is given in the following formulas (1) and (2). 

I '■• io= r • F - u ' . ^ 

m ♦ x = — k * x - c ♦ x — F (2) 

wherein I is the wheel inertia! moment (e.g., 0.8 Kgm 2 ), m is the mass of the 
wheel (e.g., 10 Kg), r is the effective radius of the tire (e.g., 0.3 m), kis the spring 
constant (e.g., 2.0 x 10 s N/m), c Is the damper constant (e.g., 100 Ns/m), F is the tire 
' friction force;, u is the braking torque inputinto the suspension - tire resonance model, 
cd is the wheel angular velocity, and x is the unsprung longitudinal displacement 

When the tire friction force F is given by the following formula (3), the transfer, 
function Gwu (s) from the braking torque to the wheel angular velocity is given by the 
following formula (4). 

# * 

x — r • oi 
F = k « • N • r^- 



(3) 



wherein kt is the u gradient; N is the wheel load (e.g., 3000 N), and V 0 is the 



vehicle speed (eg., 10 m/s). 

r . s "v. -s' + C Nk, +cv, ) - s -f-kv. 

Cs)=— ■ 

■Iv. - s'+CNk, I+cv.Ifnr'Nk, )• s'Kkv.I+cr'Nk . ) ♦ s+ kr J Nk« (4) 
Rgs. 3A and 3B are Bode diagrams of formula (4) when the u gradient kt is 2, 5, 
10, and 20. It can be understood from Rgs. 3A and 3B that as the ll gradient 
increases and a margin is generated in the tire characteristics, the characteristics begin 
to show oscillation. 

Next, the u gradient estimation will be explained, the ll gradient estimation 
may be carried out on the basis of formula (4), however, this explanation will use the 
second order resonance model of formula (5) (a simplification of formula (4)) which 
makes the unsprung resonance frequency a natural frequency. 
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wherein K is normal gain, con is the undamped natural frequency, and C is the 
damping ratio. Further, instead of directly performing an estimation of the ll gradient, 
the damping ratio ^ of this model is identified and the ll gradient is estimated from 
the identified damping ratio £ . The il gradient k t and the damping ratio Z, form the 
approximate relationship given below which represents an inverse number. 
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Accordingly, by identifying the damping ratio, it is possible to estimate the u 
gradient using the above forrhula.The resonance characteristics in the natural 
frequency are clearer the smaller the damping ratio of the second order resonance 
model, therefore, when the damping ratio is identified as being small, the n gradient 
is estimated to be large. Conversely, when the damping ratio is identified as being 
large, the n gradient is estimated to be small. . • 

The following formula (6), obtained by rearranging the wheel velocity and 
braking torque from the above formula (5), is used to explain the identification of the , 
damping ratio. 

2 C <o a sw— — s*w — o>„ 2 w+K <o n 2 u (6) 
In formula (6), as is explained below, the braking torque u, which is an input 
factor, can be ignored, the formula (7), which is obtained by ignoring the braking 
" torque, is made discrete, and the damping ratio can be identified by applying an online 
identification method to the formula (11), which is obtained by making formula (7) 
discrete. 

Firstly, if the input of the braking torque is ignored, formula (6) is expressed by 
formula (7) below. 

2C o>„ sw=^s 2 w -a>„ 2 w (7) 
10 



Formula (7) is made discrete by substituting the Justin conversion of formula (8) 
below in formula (7). 



2 ( 2 - 1 ) 

s = 



r (z + 1 ) (8) 
wherein r is the sampling time. 

Note that the Euler's backward method expressed in formula (9) below may be 
used as the method for making a formula discrete as well as the "Justin conversion 
method. 

z - 1 

S = ; 

T Z 

. ■ - . . - O) 

Formula (7) resulting from the Tustin conversion can be rearranged as in formula 
(10) below. 

4eo» t (2; 3 -l)wC=(o) a 2 t 2 +4)(z-*+1)w+(2 <o n 2 r 2 -8)w ( 10) 
Formula (10) can be expressed as a time domain giving formula (11) below. 
0 C = y 

wherein 

0=4 o> n r 0r[i-2j-»r[i]) 

y=( o>„ 2 r 2 +4)(w[i-2]+wCiJ) + (2< U „ 2 r 2 -8) w [i-lj 

(11) 

Here, formula (11) has been obtained by making formula (7) discrete and, 
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therefore, formula (11) corresponds to formula (7). Thus the physical meanings of <t> 
and y can be interpreted as follows, <J> is the physical amount corresponding to the 
change in the wheel velocity, while y is the physical amount corresponding to the 
wheel velocity and the change in the wheel velocity. 

Next, the online identification method is applied to formula (11). Here, the 
damping ratio Z, \s identified using the instrumental variable method. The algorithm 
using the instrumental variable method is expressed as is shown in formula (12) 

below. ... 

- • z=wti-3] 

PCi-l]z 



X + tfP[i-l3z 



0[i]=0[i-l]+H - Cy— ^ * 0 [i-l3) 

PCi]=— C1-H^0 -PCi-1] 

■ ■ '■' ' X • . : .-. . • : 

; " ■ ; (12) 

wherein A is a forgetting coefficient. 

Figs. 4A to 4C show the results when the damping ratio is identified using the 
instrumental variable method. The figures show that, in accordance with the osciHation 
characteristics after the vehicle has traveled opto the dry surface, the identified value 
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of the damping ratio becomes smaller, in other words, the u- gradient becomes larger, 
and a margin is generated in the braking force characteristics. Thus, it can be seen that 
an appropriate identification has been made. 

Figs. 5A to 5C show the results of an experiment in which, at the point when the 
estimated value of the damping ratio has gone below 0.4, it is determined that the u 
gradient is large and a margin has been generated in the braking force characteristics, 
in other words, that the vehicle has traveled onto a high u (dry) road surface, and 
the braking hydraulic pressure is made to increase at a consistent rate to a pressure 
appropriate for a high ti surface (as shown by the ramp-shaped portion of the graph 
in Fig. 4B). In this way, a sufficient suitability for the road surface can be achieved. 

Note that the above explanation was for an example in which the damping ratio 
was identified using an instrumental variable method, however, the present invention 
is not limited to this and another identification method such as the least squares 
method may be used to identify the damping ratio. 

Moreover, the physical quantity estimation device of the present invention can 
also be used for determining road surface conditions or the condition of a vehicle. . 

Further, the above explanation was for an example in which the p. gradient, 
which is one of the physical quantities relating to the ease of slippage between the 
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road surface and a vehicle wheel, was estimated, however, as is also understood from 
formula (4), the transfer function Gwu (s) is also changed by variations in the wheel 
load and the. vehicle speed. Accordingly, the amount estimated may be one in which at 
least one of the variations in the wheel load and vehicle speed has been added to the 
u gradient, which is one of the physical quantities relating to the ease of slippage 
between the road surface and a vehicle wheel. 
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WHAT IS CLAIMED IS: 

1 A physical quantity estimation device comprising: 

wheel velocity detection means for detecting the velocity of a wheel and 
outputting a wheel velocity signal; and 

physical quantity estimation means for identifying a parameter of a physical 
model expressing unsprung resonance characteristics from the wheel velocity signal 
and estimating a physical quantity relating to ease of slippage between a road surface 
and a wheel from the identified parameter. 

2 The physical quantity estimation device according to daim 1, wherein the 
physical model is a resonance model which has an unsprung resonance frequency as 
an undamped natural frequency, and wherein the physical quantity estimation means 
estimates the physical quantity on the basis of a damping ratio of the resonance 
" model. 

3 The physical quantity estimation device according to daim 1, wherein the 
physical quantity relating to ease of slippage between a road surface and a wheel is a 
road surface a gradient which is a ratio between a road surface u. and a slip. 

4 The physical quantity estimation device according to daim 3, wherein the 
physical model is a resonance model which has an unsprung resonance frequency as 
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an undamped natural frequency, the parameter is a damping ratio of the resonance 
model, and the road surface m gradient and the damping ratio are in inverse relation. 

5 The physical quantity estimation device according to daim 1, wherein the 
physical model is a resonance model between a tire and a suspension in a longitudinal 
direction of a vehicle. ^ 

6 The physical quantity estimation device according to claim 5, wherein the 
resonance model is an equation expressing a relation among a physical quantity 
corresponding to a change in the wheel velocity, a physical quantity corresponding to 
the wheel velocity and the change in the wheel velocity, and a damping ratio of the 
resonance model.. ♦ 

7 The physical quantity estimation device according to daim 5, wherein the 
resonance model is a second order resonance model. - 

- 8 The physical quantity estimation device according to daim 1, further comprising a 
bandpass filter which passes on only signals of frequency bands which are related to 
unsprung resonance, wherein the physical quantity estimation means identifies the 
parameter from the signals passed through the bandpass filter. 
9 An ABS control device comprising: 

the physical quantity estimation device according to claim 1; and 



16 



control means for controlling braking pressure on the basis of a physical 
quantity estimated by the physical quantity estimation device. 

10 The ABS control device according to claim 9, wherein the physical quantity 
relating to ease of slippage between a road surface and a wheel is a road surface p. 
gradient which is a ratio between a road surface ti and a slip and the control means 
controls the braking pressure so that the road surface u gradient is 0 or . 
approximately 0. 

11 The ABS control device according to claim 9, wherein the control means increases 
the braking pressure in a ramp-shaped manner so that the braking pressure is adapted 
to a high u road surface when the control means judges that the vehicle has traveled 
onto a high u road surface on the basis of the physical quantity. 

12 An ABS control device comprising: 

the physical quantity estimation device according to daim 2; and 
control means for controlling a braking pressure on the basis of a physical 
quantity estimated by the physical quantity estimation device. 

13 A physical quantity, estimation method comprising the steps of: 

detecting a velocity of a wheel and outputting a wheel velocity signal; 
identifying a parameter of a physical model expressing unsprung resonance 
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characteristics from the wheel velocity signal; and 

estimating a physical quantity relating to ease of slippage between a road 
surface and a wheel from the identified parameter. 

14 The physical quantity estimation method according to claim 13, wherein the 
physical model is a resonance model which has an unsprung resonance frequency as 
an undamped natural frequency, and when estimation of the physical quantity is 
performed, the physical quantity is estimated on the basis of a damping ratio of the 
resonance model. • . 

15 An ABS control method in which braking pressure is controlled on the basis of a 
physical quantity estimated using the physical quantity estimation method according to 
claim 13, 

16 An ABS control method in which braking pressure is controlled on the basis of a 
- physical quantity estimated using the physical quantity estimation method according to 

claim 14. 

17 A physical quantity estimation device substantially as 
hereinbefore described with reference to the accompanying 
drawi ngs . 
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